Slow logarithmic relaxation in models with hierarchically constrained dynamics 
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A general kind of models with hierarchically constrained dynamics is shown to exhibit logarith- 
mic anomalous relaxation, similarly to a variety of complex strongly interacting materials. The 
logarithmic behavior describes most of the decay of the response function. 
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In the last years, an anomalous slow relaxation that 
can be accurately described by a logarithmic law, has 
been found in the time evolution of a variety of complex 
strongly interacting materials. This includes spin-glasses 
|],|), granular materials) , structural glasses [§-§), 
and protein models (l^Jlt]. In all the above cases, the 
decay of some perturbation <f>(t) can be fitted at interme- 
diate times by an expression of the form 



<f>(t) = d-C 2 lnt, 



(1) 



where C\ and C2 are constant, and may depend on the 
macroscopic variables (temperature, compactivity, . . . ) 
characterizing the state of the system. The time range 
in which the logarithmic law accurately describes the re- 
laxation covers most of the change of <fr(t) from its initial 
value to the final one. 

Palmer, Stein, Abrahams, and Anderson (PSAA) in a 
seminal paper Jt2| described some dynamical models for 
relaxation in complex systems with strong interactions. 
The models were based in a hierarchically constrained 
dynamics so that the relaxation of the system involves a 
sequential series of correlated processes. Although PSAA 
concentrate in those physical systems showing "stretched 
exponential" relaxation, it is clear that the general pic- 
ture emerging from hierarchically constrained dynamics 
can be translated into very different physical situations. 
As a consequence, the above models provide a very useful 
tool to investigate many kinds of "anomalous relaxation" 
in complex systems. 

In particular, the picture of hierarchical constrained 
dynamics seem to be physically appropriate to describe 
many of the complex systems showing logarithmic re- 
laxation. Let us consider as a particular example the 
densification of powders and structural glasses at high 
hydrostatic pressure studied in Refs. ||,|]]7]||. External 
pressure causes some selected mesoscopic regions of the 
system to reach a more compact state very quickly. The 
nearest neighbor regions are influenced by this rearrange- 
ment, and can also reorganize themselves. Of course, 
these new structural changes will affect the regions out- 
side the nearest neighbors, allowing them to relax, and 
so on. The above discussion suggests that the relaxation 
in this kind of systems is slow because it consists of a 
large number of correlated reorganization steps with an 
increasing characteristic time. In other words, it has a 



hierarchical structure, with the faster degrees of freedom 
(associated to the "inner" regions) controlling the relax- 
ation of the slower degrees of freedom (associated to the 
"outer" regions). A similar picture has been already sug- 
gested in Refs. |||| . The aim of this paper is to introduce 
a quite general class of hierarchical models, which seem 
appropriate to describe this kind of systems, in order to 
derive the logarithmic relaxation of Eq. (Q). 

Like PSAA we will formulate our model in the most 
general way, hoping it characterizes the behavior of a 
specific class of systems. We assume that the state of the 
system can be given by means of N modes. They do not 
correspond to a particle description of the system, but 
to some mesoscopic level of description. In this sense, 
N does not scale with the size of the system, but it is 
assumed to have a definite value, independent of the size, 
for a large enough (macroscopic) system. Of course, the 
specific value depends on the system and the physical 
situation we are dealing with. 

The N modes can be grouped into a discrete series of 
levels, n = 0, 1, 2, • • ■. The number of modes belonging 
to level n will be denoted by N n , so that N = J2 n ^n- 
Moreover, modes in level n will be indicated by the vari- 

(n) 

ables Si , i = 1, 2, • • • , N n . For the densification exper- 
iments of Refs. @,||,@,D, we can regard S^ as a pseu- 
dospin variable, with its two possible values correspond- 
ing to different volumes of the mesoscopic regions. The 
dynamics of the system is defined in a hierarchical way as 
follows. A mode in level n + 1 can only modify its state 
if the configuration of modes in level n belongs to a well 
defined subset of all the possible configurations of the 
modes in that level. Clearly, this slows down the modes 
in a given level, as compared with those of the previous 
level. More specifically, the average transition rate Wn±i 
for a mode in level n + 1 is assumed to have the form |l2] ] 



W n+1 = W n P n . 



(2) 



where P n is the probability that the modes in level n are 
in any of the configurations allowing to relax level n + 1 . 
Our choice for P n will be 



P n = cxp(-JV„A/i), 



(3) 



with A/i being a dimcnsionless free energy defining the 
activation barrier per mode. Eq. (||) appears quite natu- 
rally if P n is evaluated from the equilibrium distribution 



1 



of the system under consideration, provided that the sub- 
set of configurations of level n allowing to relax level n+ 1 
involve all the modes (or an extensive part of them) in 
level n. Of course, the free energy A/i must be appro- 
priately interpreted in each particular situation. For in- 
stance, in the densification experiments at high pressure 
11110]!] it would be an activation enthalpy divided by 
the temperature of the sample. On the other hand, in a 
granular system, it could be understood as an "effective 
volume" divided by the compactivity |p^|-|l6||. Thus, the 
corresponding average relaxation time of the mode is 



'n+l 



,exp(iV n A/i). 



(4) 



For the sake of simplicity, the barrier Afi is assumed to 
be independent of the level n of the mode and also of the 
number of modes in the level N n . Iteration of Eq. (Q) 
yields 



exp ( N„ 



(5) 



The range of variation of n and the values of the pop- 
ulations N n characterize the complexity of the time evo- 
lution of the system. In principle, there is no need for 
a maximum value n max of n, i.e. n max may be infinite. 
On the other hand, the systems we are interested in show 
slow relaxation, which requires a large relaxation time. 
In this way, it is necessary that 



lim T n > Tq . 



(6) 



This condition for slow relaxation was already pointed 
out by PSAA @. 

We are interested in the relaxation of the quantity 



Q(*) = ^££$ B) (*)>. 



(7) 



where the angular brackets denote configuration average. 
Note that in the context of the densification experiments 
H[l@]i] Q(t) would be proportional to the actual volume 
of the sample at time t, with the physical meaning given 

(n) 

above to the variables S { . The corresponding relax- 
ation function is defined in terms of the asymptotic long 
time value Q(oo) as 



Q(f)-Q(oo) 
Q(0)-Q(oo) 



^ w„ exp 



(8) 



with io n — N n /N being the fraction of modes in level n. 
Combination of Eqs. (ph and M) gives 



n— n max 

<j){t) = uj n exp 

n=0 



n-1 



— exp -C V l 



(9) 



where 



(10) 



Next we will assume that u n changes very smoothly with 
n. This allows to introduce a continuum limit for the 
distribution of relaxation times, which is expected to 
be closer to real systems than the discrete level model 
considered up to now. We define a function uj(x) by 
cD(ne) = cj„, with e < 1, so that the variable x n = ne 
is almost continuum and the sums over n can be re- 
placed by integrals. Let us introduce the new variable 



dx'ujx 1 ) 
J Q Xmax dx'Q(x') ' 



(11) 



< u < 1, which measures the fraction of the total 
number of modes N which are in levels up to a value 
n = x/e. It is then a simple task to transform Eq. (|p) 
into 



du exp 

(i V TO 



t 



(12) 



showing that the decay of cj> can be characterized by only 
two parameters, namely tq and £. The former fixes the 
time scale, while the latter determines the actual shape of 
the relaxation. The elementary relaxation times given by 
Eq. (||) are t(u) = To exp(£it), and, in particular, r max = 
To exp £. Thus, the condition for slow relaxation, Eq. (|^), 
implies that exp ( > 1. 

The analysis of Eq. @ is greatly simplified by realiz- 
ing that it is equivalent to 



TO 



tq 



(13) 



where E\(z) is the integral exponential function [ jl7| . 
There are two natural time scales in Eq. (p^|). One is 
that defined by t/ro, while the other is a slow time scale 
defined by s — — — . In the following, we will 

use the asymptotic expansions of the exponential inte- 
gral function from Ref. ]l7| ]. For very short times such 
that t/r < 1, Eq. (ph is equal to 



1 



< t 



C to 



O 



(14) 



i.e., (f>(t) decays very little from its initial unity value in 
this region. In the opposite limit of very large times for 
which s 3> 1, it is 



1 e~ l 



(15) 



the response function has already decayed to very small 
values. One main conclusion of the above discussion is 
that the relevant part of the relaxation, in which <p(t) 
changes from values very close to unity to very close 
to zero, takes place between the two asymptotic time 
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regimes we have analyzed. In this way, we are led to to 
consider the distinguished limit t/r > 1, s < 1. In this 
intermediate time window it is found that 



<t>(t) ~ 1 - j ( 7 + In - 



(16) 



4> 




FIG. 1. Relaxation function (f> plotted as a function of 
time. From left to right, the curves correspond to exp £ = 10 3 
and 10 5 . The integral representation of <j>, Eq. (|l^), is given 
by the circles, while the solid lines are the asymptotic expres- 
sion for intermediate times, Eq. jl^). 

This is precisely the logarithmic decay of Eq. (jl]). 
Moreover, the previous discussion of the short and large 
time limits, clearly indicates that most of relaxation of 
the system is expected to be accurately described by Eq. 
(|l6|). As an example, we present in Fig. [l]the relaxation 
function <f>{t) for exp£ = 10 3 and exp£ = 10 5 . The inte- 
gral expression given by Eq. (jlj) (circles Hs well htted by 
the logarithmic asymptotic form of Eq. (16) (solid lines). 
The agreement extends over several decades of interme- 
diates times, measured in units of the basic time tq. The 
fit improves as the value of £ increases, consistently with 
the condition for time scales separation. 

It is also interesting to compute an average relaxation 
time f, given by 



dt (j)(t) = To 



£-1 

c 



(17) 



showing an exponential increase with £. Taking into ac- 
count the definition of (, Eq. (10), the above result can 
be understood as corresponding to an activation energy 
barrier whose height is proportional to the complexity 
of the system, measured by the total number of modes 
N actively involved in the relaxation. As already men- 
tioned, for a given specific situation, the dimensionless 
parameter Afi might have a well defined physical inter- 
pretation, being for instance the ratio between the Gibbs 
free energy barrier per mode and the thermal energy. In 
this way, it is also interesting to note that the coefficient 



of lnt in Eq. Jl6| ) is C _1 j which is then proportional to 
the temperature T. This would lead to the Tint behav- 
ior observed in many physical systems, for instance in 
Refs. §§. 

In summary, we have studied a simple class of models 
which exhibit an anomalous logarithmic relaxation. The 
dynamics is hierarchically constrained, and they are ad- 
equate to study slow relaxation because there is a clear 
separation of time scales. The slow logarithmic relax- 
ation is similar to what is experimentally observed in a 
variety of complex systems, fitting the data in most of 
the relevant part of the relaxation, i.e., the time window 
where the relaxation function is neither too small nor too 
close to its initial value. Perhaps the main assumption 
of this work is the "extensiveness" with N n of the effec- 
tive barrier for modes in level n + 1, as expressed by Eq. 
(||). This is quite a general condition and, therefore, the 
class of hierarchical models presented here may provide 
an explanation of the appearance of linear logarithmic 
relaxation in very different physical situations. In this 
context, the role of hierarchically constrained dynamics 
may also be relevant in order to understand other kinds of 
anomalous relaxation, as the inverse logarithmic behav- 
ior found in compaction experiments of granular systems 
p8| , ^9| . Also, it seems interesting to investigate the be- 
havior of hierarchical models in the temperature cycling 
experiments usually carried out in glassy systems ]2C|]. 
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